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ABSTRACT 


Backscatter of acoustic signals generated by surfaces 
having point scattering and line scattering features is 
studies experimentally in air. Signals are analyzed for 
propagation over three models: a smooth plane wedge, a 
wedge whose surface is slightly roughened with uniformly- 
packed lead shot, and one randomly roughened with gravel. 
Scattering from the rough surfaces at normal incidence is 
measured to estimate the statistical parameters that describe 
the surface roughness. The returned acoustic energy from 
the composite area of rough surface and wedge for low grazing 
angles is analyzed to determine whether the identity of the 
two individual backscatter effects can be separated in their 
combined spectrum. The near grazing backscattered amplitudes 
are much larger than predicted by conventional rough surface 
Scattering theory and suggest the presence of a backscattered 


boundary wave. 
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fie  LNTROBUCTION 


Diffraction of acoustic energy from both smooth wedges 
and rough plane surfaces has been theoretically investigated 
by a number of authors. Many at sea measurements have been 
reported for rough surfaces. Of particular interest is the 
theoretical development of Biot and Tolstoy [1] for impulse 
diffraction of point source radiation by plane rigid wedges, 
and the experimental results for sea surface scattering, 
eeg., Clay and Medwin [2] and Urick [3]. 

Biot and Tolstoy theory has been tested experimentally 
by Bremhorst [4] and Spaulding [5] for simple wedges and 
complex seamount modeling, respectively. The experimental 
results showed close agreement with theory. 

This study continues the experimental comparison of 
diffraction by various effects in the specific case of 
backscatter. The approach involved determining the spectral 
relationship for the smooth wedge and a rigid, random rough 
surface, separately. A wedge covered with a rough surface 
of gravel was used to attempt to differentiate between the 
returned signals from both the wedge and rough surface 
effects. Also, a wedge covered with lead shot was used to 
present a uniformly-packed, slightly rough surface. The 
results were intended to show how the combination of wedge 
slope angle and surface roughness interact. This relationship 


would ultimately identify a range of these characteristics 





where one Signal is not extremely dominant over the other, 
but rather the two signals could be deconvolved to show their 
separate effects; and thereby identifying the character of 


the source of backscatter. 





Mee kcoomnRocn eA CLL ETIES 


A. ANECHOIC CHAMBER 

The experiments were conducted in the Naval Postgraduate 
School anechoic chamber (air), which provides excellent iso- 
lation from external noise sources, absorption of internal 
sounds, and atmospheric stability. The internal design of 
Fiberglas wedges allows for absorption of approximately 993 


of the incident sound of frequencies greater than 100 Hz. 


B. DATA PROCESSING EQUIFMENT 
Data acquisition and processing were accomplished by 
using a digital computer system composed of four major com- 
ponents interfaced to provide high speed analog to digital 
conversion, digital processing, and data printout. The 
design was developed by the Special Projects section of 
Naval Air Development Center in conjunction with Pinkerton 
Computer Consultants, Inc., of Warminster, Pennsylvania. 
The major components are: 
im Interdata Model /0 Computer 
This minicomputer is a digital design that is FORTRAN 
and BASIC programmable with a 64 thousand byte memory. Addi- 
tionally, data stored on floppy discs or digital tape 
cassettes can be input for processing. 
eee noenix Analog to Digital Converters, Modell ADE 912 
iwoeA De. 912 analog to digital converters may be used 


separately or simultaneously. Each converter is a high 








speed, high accuracy device capable of encoding +10 volt 
input signals into twelve binary bits of data with a reso- 
lution limit of one part in 4,095 at a maximum rate of 2 
microseconds per conversion. It measures the input voltage 
relative to an internal precision reference voltage with 
pummerecuracy Of +0.025% of full range. The ADC 912's will 
accommodate a typical commutating through-put rate of 
476,190 channels per second, including settling time. The 
Sampling frequency is set ina stable frequency oscillator 
and sent to the converters via a sampling circuit. As used 
conservatively at NPS, the maximum sampling frequency is 
bzenvo0 Hz. 
3. Orbis Diskette Drive, Model 76S 

The Model 76S diskette drive is a peripheral device 
interfaced with the Interdata 70 Minicomputer in order to 
read and write on an industry standard floppy disc. The dual 
drive unit provides 3.2 million bits of data with a data rate 
of 250,000 bits per second encoding for each drive and is 
capable of a 6 millisecond access time track-to-track. 


4. Texas Instruments Silent Electronic Data Terminal, 
Model 733 


The TI 733 consists of a keyboard used as a program- 
meng input/output control device, a printer, and transmit/ 
receive mechanism for the peripheral disc drive unit. The 
overall system facilitates rapid, accurate processing of any 
desired type of analog electrical signal and was used pri-~ 
marily for frequency domain analyses using standard FFT 


algorithms. 
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C. STANDARD EQUIPMENT LIST 
The following scientific equipment was utilized to 

conduct the experimental work: 

Interface Technology Timing Simulator/Word Generator, 
Model RS-648 

Wavetek Arbitrary Waveform Generator, Model 175 

Instruments, Inc., Power Amplifier, Model L2 

E&L Instruments DD-1 Digi Designer 

Tektronix Type 551 Oscilloscope with two dual-trace 
plug-in units 

General Radio Type 1312 Coherent Decade Oscillator 

Hewlett-Packard Electronic Counter, Model 5223L 

Fluke Digital Multimeter, Model 8000A 

Hewlett-Packard DC Power Supply, Model 721A (two) 

Princeton Applied Research Pre-Amplifier, Model 113 (two) 

Krohn-Hite Frequency Filter, Model 3550 

Bruel and Kjaer Type 4145 Condenser Microphone 

Bruel and Kjaer Type 4133 Condenser Microphone 

Bruel and Kjaer Microphone Preamplifier, Model 2619 

Bruel and Kjaer Microphone Power Supply, Model 2804 

The equipment diagram for the experiment is shown in 


proure 3. 
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Figure 2. Typical Rack Mounting 
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a On 


Scattering can be generally explained as the phenomenon 
which occurs whenever a wave front encounters an obstacle, 
thus producing a new source by the principle of Huygens wave- 
lets. With the exception of forward scatter, we can separate 
the sound waves scattered by an object from the direct sound 
wave by uSing a pinging technique and observing the scattered 
Signal when there is no signal from the source [2]. This 
experimental work involves the study of two theories for 
backscatter--that for which sound impinges on a change of 
slope such as for a wedge and that for sound scatter from a 


moe Surface. 


A. WEDGE 

To examine the former theory, consider the model of an 
infinite wedge which is bounded by rigid plates at 6 = 0, 
6 = oe Me nteeeectiing along the apex, chosen as the z-axis. 
The region is occupied by a homogeneous compressible fluid 
of sound speed, c, and density, op. 

To determine the normal modes of the system, it is 
convenient to solve the acoustic wave equation in cylindri- 


cal coordinates, in terms of the displacement potential ¢: 


a pi eee ait 3.5 
sre or ee 397 - ae 5 


a 
—- 
are 
a 
Ss 
i 
“—- 


) 


Kite 


. 
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The harmonic solutions to this equation are of the form 


en a. nae 


Sy? oa 
(r,0,z,t) ape (2) 


(xr)e 
where v, «, and y are separation constants and may be thought 
of as components of the wave propagation vector, k, in the 


azimuthal, radial, and z directions, respectively, and 


Boundary conditions can be established such that the 
Sides of the wedge are assumed to be perfectly rigid, re- 
quiring the displacements normal to the sides to vanish, 


and can be expressed mathematically as, 


3b _ Z 
a Gee at o—200eo QO... (4) 


Writing the first term of the RHS of equation (2) as 


+ivée 
e = cos v8 + s1n v8 


and applying the boundary conditions yields the solution 


Bemegquation (1), 


(12) etiyzZ ativt 
Vv 


(5) 


6 = cos v_8H 
n 


eG 





where 





v = —. (6) 


At r = 0, the imaginary part of the Hankel function 


becomes infinite, which is not physically realizable. 
E20 


Therefore, only the real parts of H 


, which are Bessel 
functions of the first kind, will be kept. 
Biot and Tolstoy [1] develop this theory by assuming an 


explosive point source at r and z = 0, which leads to 


0’ °o! 
solutions to equation (5) that are symmetric in the z-axis. 
Applying the method of normal coordinates, they arrive at 


the solution, 


98 | dt 
ae om L cos v @ cos v6) I, , ct >z (7) 
where 
T= fd) (er)T, (erg) Tle lett = 2°) 41 ede (8) 
0 n n 


For a detailed look at the mathematics of the solution tech- 
nigue, see Refs. 1 and 6. A transform given in Appendix X 
of Ref. 6 may be used to reduce I. to more physically com- 
prehensible terms. 

Now let 


Li 2 (9) 


ct 
ll 


ip 2 Z 
al(r-25p) +22, ] 


Le 





a ti 7 





t' = =[(r+xr,)* +2774 


where to is the least possible time required for a sound 
pulse to travel directly from a source to receiver and t' 
is the time required to travel from source to wedge apex to 
receiver. Equation (8) will have three solutions: t < tos 
i 


foe Gl, hvG!lUcS cE. [6CFOr El U< ct I. = 0 as would be expected 


0 0’ 
and there 1S no signal. For to Selo eG | ene se uewon des— 
cribes the acoustic field for the direct energy (not of 
interest), reflected energy, and any backscattered energy 
prior to the signal from the wedge apex (for the rough sur- 
mteemedse Considered). For t' < t, the solution is the one 


of interest for the theory applied to backscatter from wedges. 


For this case, 


ale ~Y ny 

ie = = SIN VT CUM (10) 

n mr Lp Sinh y n 
where 

o*t* = (x* trp* + 2°) 
y = arc cosh ae Cie) 
0) 

USing this solution, equation (7) becomes 
a= AS Bees cos v_6@,cosv_6 sin v re ™ , (12) 
dt 75,0 rr,sinh y in 0 n n 


= et 
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Rewriting the trigonometric functions in terms of their 
exponential identities and then collecting the conjugate 
pairs, equation (12) is rewritten as 


1 TY / ey 


¥ reXo Sinh y 7 (13) 


a 
o 
Q 





@ 
(+ 
i 
=| 

@D 


Sin(1/8,,) (126269) 


“TY/ 81, ~27y/6,, 
1-2e cos (7/6,,) (n+6+68,) +e 


} 


when the + signs indicate that the term in brackets is actualy 
a sum of four terms corresponding to all possible combinations 
of these signs. 

For the unit impulse source used in Refs. 1 and 6, the 
displacement potential 4 is related to the acoustic pressure 


by 


Q 
oS 


. + eS (14) 


Q> 
(t+ 


If, on the other hand, one assumes a source that is a delta 
function in time as well as space, the acoustic pressure 


resulting from the backscattered wave alone is given by 


Mme ec te 
an Ory rry Sinh y 
Sin(7/6,)) (26269) , os 
a) ae ~2ny/8,, 
i1-2e cos (1/6,)) (n2646,) +e 


ILS 











where 


oft* = (ro +r * 42%) 
y = are cosh 
ee is 
0) 
and 
a. OF fluid filled region above the wedge 
8) = angle between the wedge side and the source 
6 = angle between the wedge side and the receiver 
ry = distance from source to wedge crest 
te = distance from wedge crest to receiver 
z, = distance along wedge crest measured from 
point of intersection of least time travel 
path 


See figure 4 for illustration of geometry. 

The backscattering process can be described as follows. 
A point source at a distance Lo from the wedge crest trans- 
mits a spherical wavefront. The wavefront first impinges 
the crest along the least time travel path at time, t = Lo/c- 
The wavefront then continues to arrive at progressively later 
t along increasing distance z. This in effect produces a 
continuous series of point sources along the crest which re- 
radiate. Thus the crest of the wedge acts as a line source 
with time shading along the line. 


Since the frequency domain solution to the impulse wave 


Scattering problem is more familiar to the applied acoustician, 
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Figure 4. Typical wedge geometry 


as 


equation (15) has been transformed to the frequency domain 
by Medwin [7], and the computed theoretical values have been 


compared to experimental data. 


B. ROUGH SURFACE SCATTER 

The statistical approach to the scattering theory des- 
cribes the random surface by the distribution of its sur- 
face heights about a mean level and by a correlation function 
that relates the heights at various points on the surface. 
The distribution of surface heights is assumed to be Gaussian 
with zero mean and standard deviation, o, which is identically 
the rms surface height. The correlation function is assumed 


to be 
2 
C(t) = exp <-(t/L) ( 156)) 


where L is the correlation distance of the surface and +t is 
the distance between any two points whose correlation is 
desired. 

The derivation of this theory is given by Beckmann and 
SpizzZichino [8] for electromagnetic radiation. However, 
acoustic scattering is analagous to electromagnetic, where 
the pressure field is developed by using the Helmholtz inte- 
gral and applying boundary conditions in conjunction with 
the Kirchhoff approximation [2]. The further assumption 
that the incident pressure is constant over a rectangular 


insonified area and zero elsewhere is imposed to yield 
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eeceormula for the scattered pressure from a random rough 


surface, 


_ F iver 
Eee hogs ac dxdy bys) 
Pq = specular scatter from smooth surface 


ik Pig AGos a4 


P59 = rr | exp (ikr,) (18) 
Pip = incident pressure 
54 = angle of incidence 
85 = angle of scatter to receiver 
ro = range from surface to receiver 
- 1 + cos Oy cos 9, - sin ay sin 85 cos 4 
ae cos 8,(cos §, +cos 8.) 7 
1 dl 2 
Vv = propagation wave vector, a function of 


geometry. 


Since P5 is a complex quantity, the mean square scattered 


pressure <p *> must be calculated in order to introduce 


2P2 
the scattered intensity. This then is related to the mean 


and variance of the scattered pressure by the expression 


<P5P5*> = <P5>°<Py*> + Varlp,] : (20) 
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The mean scattered pressure is equivalent to the coherent 
component of the specularly scattered pressure field since 
it has constant phase with respect to the smooth reflected 


pressure. 


2° = Pog F exp(-g/2) sinc vx sinc hi (28) 
where 
Sin v X 
X 
ane. SS 
x vi 


and g is defined as the roughness parameter, 


g = k 707 (cos 6, + cos 8) ° ‘ (22) 


The variance is derived by assuming surface roughness to 


Bemisotropic, 


T Fe L° exp (-g) oa 
<> 
= gm - (v7 v4) EB 
L mm! eae 4m : ee 


and this equates to the incoherent component of the scattered 
pressure field. A more detailed analysis of the mathematics 


in this section can be found in Ref. 8. 
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For the specular case, the scattered intensity reduces 
to 
m 


3] (24) 
, mm! 


<p p.*> 2 
ot = exp(-g) + T= [exp (-g) 


ee 


ll o~7 8 


mM 


The expression in brackets is defined as S(g), and approaches 
a constant linear slope at very low or very high values of 


roughness, 


Sg). —C fois ) Geo. 1 
Sy aly Gy Ot ec eo ak 


For normal incidence, the relative scattered intensity 
Shows a relationship to frequency squared through e J for 
g < 1 and approaches an asymptotic value of Payee Ou 
eee «LO. 


In the backscatter case, a Surface backscattering factor 


is defined as 


<p.p x e 
2 pee ob. es (25) 


S 
B 2 
1P1 | 


where p, is the incident pressure at the surface. love el calle e 


mevely smooth surface, g << 1, the intensity becomes 


pine 
US am te ey 


<PoP5*? = $48 — exp (-k*L*sin* 89,79) + ee) 
eneCOS ©2 


2 





and the backscattering factor is 


Be 2 
= aot exp (-k“L- sin“9, -g) a2 
A 


For the rough surface, g >> 1, the coherent term 1s considered 


negligible so that the variance equals the intensity 





S Zane 
< * > Varl[p.] me P20 © xD ila tty Saag ] 
— = ee ee SS = 
P4P> Po Ag p y g 
(28) 
and the backscattering factor becomes 
= ~ Vee cos*@ van exp [-tan- 9 Vine (29) 
B at 1 
where 
2 oe 
E = s = mean square Slope . (3.0) 
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IV. EXPERIMENTAL DESIGN AND PROCEDURES 


Pema lLGN OF THE MODELS 

The choice of models to experimentally test the theory was 
a progression from a simple wedge to more complex roughened 
surfaces and finally to a combination of these two features. 
It was necessary that the models appear to the source and 
recéiver as being spatially infinite compared to the wave- 
length of the signals. By using sufficiently large model 
dimensions and a pulsed signal technique, interference from 
edge diffraction was avoided and the surface could be con- 
feegerea infinite for the frequency band of the experiment. 
Additionally, the models were required to approximate a per- 
fectly rigid condition. Aluminum was selected for the build- 
ing material due to its availability, machinability, rigidity, 
and reflection coefficient in air of almost unity. 

f. Smooth Wedge 

The wedge fabricated for Spaulding's [5] experimental 

work was available and appropriate for this part of the study. 
The dimensions--two pieces of 1/4" thick aluminum, each 1.52 
feemeemert.) long by 0.6 m. (2 ft.) wide, joined at the crest 
BPemming an interior angle of 152°--were dictated by the size 
of an aluminum section available, which satisfied the infinite 
surface criterion. The symmetrical angle of 14° with the 
horizontal was chosen to approximate the average upslope 


and downslope of the seamount modeled by Spaulding. For 
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this study of backscatter, the wedge presented a depression 
angle of 28° as the sound wave traveled up the slope and 
pmeauntcered the effect from the crest. 
Z. Rough Surface 

The rough wedge deSign required a random surface 
whose rms height and rms slope of roughness were known 
parameters for comparison to theory. Since the final model 
was to be a combination of a wedge and rough surface, the 
rough surface was prepared on a model of identical dimensions 
as those for the smooth wedge. Due to weight and cost, 
aluminum of 3/16" thickness was chosen, yet the massiveness 
was maintained by the addition of lead shot as the roughening 
element. Since the wedge had a dual use with another experi- 
mental study, lead shot was chosen to be compatible for both 
experiments and was distributed uniformly in a dense packing 
pattern over the surface. It was intended that the frequency 
range of the study be kh < 1 where h is the separation of 
the shot so that the regularity of the uniform packing would 
not affect the backscatter. 

However, for the wavelengths available (due to the 
Kh < 1 limitation), the size of the lead shot and the uniform 
packing presented a surface that appeared "acoustically 
mimeen 1.6€., g < 1. Furthermore, since the theory for a 
random rough surface assumed GausSian distribution for the 
roughness, a different roughening element was required to 


approximate this distribution. The wedge-shaped base remained 
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the same with gravel applied to roughen the surface. The 
Statistical quantities of this surface were determined from 


height measurements by Lt. E. Jordan. 


B. SOURCE/RECEIVER SELECTION 

Bremhorst [4] conducted extensive experimental testing 
to determine the most suitable source and receiver to satisfy 
the requirements of size, frequency response, acoustic out- 
put, and receiver sensitivity. After a compromise of these 
parameters for the best experimental combination, Bruel and 
Kjaer (B&K) condenser microphones were chosen as both source 
and receiver. 

To examine the Biot-Tolstoy theory, it 1s desirable to 
approximate a point source to transmit spherically symmetri- 
cal waves. This condition requires that the source dimen- 
Sions be smaller than the wavelength of the transmitted 
Signal (ka << 1, where k is the wave number and a is the 
radius of the source). To approach the assumption of a point 
source, the smallest possible microphone should be used as 
a source; however, the limiting factor for size is the need 
for a usable acoustic pressure output over the frequency 
range of interest. This was especially critical when trans- 
Mitting over the rough surface, as the effects of secondary 
scatter added to the transmission loss. Although the B&K 
condenser microphones operate as piston radiators, they 
follow the fundamental property of radiating a spherical 


diverging wave but with a directivity factor included. The 
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directivity term could be ignored as long as the sound 
field was essentially constant over the insonified area of 
interest, i.e., on or near the acoustic axis of the 
microphone. 

The receiver selection involved the trade-off between 
the smallest possible microphone for good high frequency 
response and the largest receiver to maintain a sensitivity 
level sufficient to detect the weak, returned signals. The 
optimum combination was using the 1" B&K 4145 as the source 
and the 1/2" B&K 4133 as the receiver. 

Another important consideration was that the pulse length 
of the single-ping signal determined the resolution of the 
Surface roughness that could be recognized. The B&K conden- 
ser microphone exhibited the disadvantage of transforming 
a Single-ping waveform into a repetitive periodic waveform 
decaying with time, or a "ringing" effect. The result was 
a larger area of insonification than desired because of the 
increased time of the output Signal. 

This study embarked on the task of constructing and test- 
ing an alternate microphone that could generate the required 
acoustic power while maintaining receiver sensitivity, yet 
reduce or remove any ringing effect. A solid dielectric 
Mylar transducer was designed using the guidelines of Ref. 9. 
Four transducers were fabricated--two 3" and two 1", with 
Peemesize having different degrees of roughness of the back- 


plate. It was projected that if these designs decreased the 


6S 


ringing effect, a newer material of polyvinylidene fluoride 
(PVF. ) could replace the Mylar surface to further enhance 

the output capability and the receiver sensitivity. This 
material has the advantage of being prepolarized to allow 

more electrical signal to either generate or detect the acous- 
elemsignal. 

Testing indicated that the Mylar and PVF, dielectric 
transducers also showed a ringing effect, so that three 
additional tests were conducted to compare the amplitudes 
and duration of received signal of all the transducers. A 
common source was used to radiate to the various receivers, 
and likewise a common receiver was used to compare the 
transducers as sources. The sources and receivers were 
facing each other on their respective acoustic axes. The 
third test used a common receiver with the various sources, 
but compared the reflected signal due to normal incidence 
on a rigid plate (Table l). 

A 3" Mylar transducer was used as a common source with 
PMiemiG@esulting differences in receiver characteristics being 
greater than any of the other tests. The 1" B&K microphone 
Meeuree) had a signal strength 23-29 dB greater than any of 
the dielectric transducers, but demonstrated a ringing of 
600 microseconds that was from two to four times longer than 
the others. - 

With an 1/8" B&K microphone as a common receiver, the 


sources were more similar in their results. The 3" Mylar 
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Transducer Comparison--All values determined 
from oscilloscope photos 


A. 3" Mylar source 


RCVR AMPLITUDE [dB re lv] DURATION [usec] 
1/8" B&K -18.4 150 
1/2" B&K - 1.4 300 
1" B&K 0.8 600 
mae), -22.0 300 
1" Mylar =26.0 200 
Saeetivlar -24.4 150 


B. 1/8" B&K receiver 


SOURCE AMPLITUDE [dB re lv] DURATION [usec] 
iy2° BE&K Salis 50 
1" B&K -14.0 100 
ay PVF, =L6<5 50 
iemviylar -14.0 40 
3" Mylar =O 25 


C. Reflection off smooth, rigid plate at normal incidence; 
1/2" B&K receiver 


SOURCE AMPLITUDE [dB re lv] DURATION [usec] 
i B&k - 3.1 180 
i” PVF, -18.8 160 
ieeMy lar -18.4 120 
3" Mylar -26.0 100 
Notes: "Duration" is measured from onset of signal to 


time when noise is greater than signal. 


a 





a i _ 


transducer had a 6 dB stronger signal than the 1/2" B&K 
and the 1" PVF, and was 4 dB stronger than the 1" B&K and the 
1" Mylar. Only a small variation in the ringing existed 
where the 3" Mylar signal duration was 25 microseconds as 
eempared to the 1" B&K of 100 microseconds. 

The reflected signal indicated that the 1" B&K source 
Signal exceeded those of the dielectrics by 15 to 23 dB, 
yet the ringing of 180 microseconds duration was only slightly 
more than the best case of 100 microseconds for the 3" Mylar. 

These results verified the known fact of poor efficiency 
for dielectric transducers. Although the B&K microphones 
were limited in the amount of applied voltage when used as 
a source, twice as much voltage applied to the dielectric 
transducers still yielded a much weaker signal. The differ- 
ence of the ringing characteristic was considered less signi- 
ficant as a variation of 100 microseconds produces a 3.45 cm. 
Spatial resolution change. Also an inherent discrepancy 
arose during the testing. The electret transducers appeared 
to produce results that decayed with repeated trials, indi- 
cating that the stability of the laboratory fabrication of 
the transducers was not on the same scale as B&K commercial 
Seoduction. 

The conclusion was to abandon use of a dielectric trans- 
ducer and revert to the previous findings--a bistatic com- 


bination of a 1" B&K source and 1/2" B&K receiver [Ref. 10]. 
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ee oLGNAL PROCESSING 
im oource Signal 

In order to obtain spectral analysis of a broad 
frequency range, the transmitted acoustic signal waveform 
must be composed of many frequencies. Ideally, an impulse 
would provide equal energy at all frequencies. However, 

a triangular waveform approximates an impulse on its posi- 
tive going pulse and provides signals at all frequencies 
with harmonics of the fundamental predominant. The harmonic 
amplitudes decay at a rate of 6 dB per octave. 

The Wavetek Model 175 Arbitrary Waveform Generator 
was used to generate the signal from stored waveforms or by 
programming any desired shape. A rise time of 500 nano- 
seconds of the triangular waveform allowed a close approxi- 
mation to the infinite slope of an impulse. The waveform 
generator was triggered by the Interface Technology Model 
RS-648 timing simulator to provide a single pulse mode. This 
method presented the ability to separate the deisred scattered 
signal from other interference for each cycle of the output 
waveform. 

The pulsed output was amplified by the Instruments, 
Inc., Model L-2 Power Amplifier to the maximum allowable 
Welecage for the B&K microphones. For the 1" B&K, a total 
AC +DC voltage of 250 V is the limit. The best combination 
of polarizing voltage and AC driving voltage was 150 V and 


100 V (peak to peak), respectively. 
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2. Received Signal 

The received acoustic energy at the microphone was 
then amplified, bandpass filtered to eliminate low frequency 
noise and prevent aliasing, and amplified again to approach 
PiemmaxlLmMum input limit of the A/D converter, +10 volts. 

The amplifiers used were PAR Model 113 preamps due to their 
very low self-noise characteristics. 

To separate the diffracted signal of interest from 
the stronger, unwanted signals of direct path travel, re- 
flected travel from the surfaces, and diffraction from the 
edges; a sampling circuit was used to select a window coin- 
ciding at the time of the diffracted analog signal. The 
Interface Technology timing Simulator accomplished the gating 
for the sampling frequency to 100 nanosecond accuracy. 

The sampling circuit was built with IC devices, asso- 
Cclated power supplies, and an oscillator. A General Radio 
Model 1312 coherent decade oscillator delivered the sampling 
frequency, since its stability characteristics provided the 
A/D converter with a constant frequency which allowed in- 
creased accuracy and repeatability of the data analysis. A 
schematic diagram of the circuit is illustrated in Figure 8. 
The circuit operates in the following manner. The decade 
oscillator porduces a constant frequency sine wave, which 
is then converted to the square wave form required by a A/D 
converter by a LM 710CN voltage comparator. The stable 


Square wave output of the voltage comparator is combined 
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With the timing simulator pulse of duration for the desired 
aperture in a SW 7408N AND gate. The resulting output of 
the AND gate is a constant sampling frequency to the A/D 
BPemvercer only for the duration of the window. The input 
pulses to the A/D converter keyed the device to sample the 
analog signal being routed to the computer through the A/D 
Seomvercer. 

Within a sampling aperture, the number of samples 
collected by the A/D converter was set to be a power of two; 
ioe O44, 128, 256, etc., as required by the FOURONE FFT 
algorithm, which performed the transform for the frequency 
domain analysis. The sampling frequency was set at an inte- 
gral multiple of the frequency of the transmitted signal. 
This relationship could be controlled by the decade oscilla- 
tor and the Wavetek 175, thereby minimizing any truncation 
of the sampled waveform and eliminating sidelobes. The fre- 
quency resolution was determined by the following formula, 


sampling frequency 


number of samples frequency resolution 


However, the limiting value of frequency resolution was fixed 
by the duration of the sampling window. 

The experiment was conducted with a 5 kHz transmitted 
Signal and a 200 microsecond sampling window, which likewise 
set a 5 kHz frequency resolution in the output data. With 
this resolution determined, then a sampling frequency of 320 


KHZ would require 64 samples of data within the sampling 
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aperture. Following one cycle of signal transmission and 
received signal processing, a 20 millisecond delay was 
incorporated in the timing sequence to allow for complete 
decay of any reverberation before the next ping was triggered. 
To improve the signal-to-noise ratio, an averaging 
routine within the computer software was performed on the 
digitized data. An ensemble of as many as 9,999 sampling 
windows could be selected for averaging prior to Fourier 
analysis. The signal-to-noise ratio could be increased by 
10 log YN, where N is the number of averaging runs. Con- 
Sidering the trade-off of computer run time versus Signal 
gain, 1,000 blocks of data were averaged to yield a 15 dB 
improvement. This averaged data was then used by the FFT 
algorithm, so that the computer output presented intensity 


Eeeespectral distribution. 
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V. EXPERIMENTAL RESULTS AND ANALYSIS 


A. BEAM PATTERN 

Prior to determining the reradiated energy from the wedge 
Or the rough surface, a beam pattern for the adjacent source- 
receiver combination was required to use as the reference 
pressure at a distance of one meter. The source-receiver 
combination (Fig. 9) was aligned to transmit at normal inci- 
dence toward a rigid, smooth plate. The source-receiver 
set-up was then rotated through oe. so that the received 
Signal described the beam pattern for the source-receiver 
combination. This procedure was performed for two arrange- 
ments that are called the symmetrical and asymmetrical cases. 
In the symmetrical case, the source and receiver are aligned 
relative to each other vertically (parallel) with respect 
to the axis of rotation, and in the asymmetrical case, they 
are horizontal (perpendicular) to the axis. As the combination 
is rotated in azimuth, both the source and receiver are the 
Same distance from the rigid reflector for the symmetrical 
Situation. However, in the asymmetrical arrangement, either 
the source or receiver is closer to the rigid plate than the 
other depending upon the direction of rotation. 

The resulting beam patterns for the two situations are 
presented in Figures 10 through 14 for various frequencies. 

Interestingly, the two patterns are quite similar regard- 


less of the relative position of the source and receiver to 
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Figure 9. Beam pattern source-receiver arrangement 
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Figure 10. Beam pattern for 10 kHz 
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Brgure il. es) Beam pattern for 30 kHz 
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Figure 12. Beam pattern for 50 kHz 
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Figure 13. Beam pattern for /0 kHz 
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Figure 14. Beam pattern for 100 kHz 
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the reflecting plate. This can be accounted for by the fact 
that the 200 microsecond aperture of the received signal 
processing 1s long enough to include the diffraction results 
from the adjacent microphone in either spatial arrangement. 
The dissimilarity between the right and left halves of the 
pattern is likewise consistent throughout the frequency 

range and results from the diffraction effects of the source- 


receiver combination and the mounting set-up. 


B. SMOOTH WEDGE 

To measure the backscatter from the crest of the wedge, 
the source-receiver acoustic axis was aimed at the crest 
(Z-axis of wedge). The position of the source-receiver 
relative to the wedge was varied to compare to the ticoae for 
the respective geometry. The position was defined by the 
Slant range from the crest to the source-receliver and by the 
angle between the acoustic axis and the wedge. The geometry 
of the experiment was limited by the interference of unwanted 
Signals when the crest backscatter could not be individually 
Gated from the other returns. Figure 15 shows an oscillo- 
scope trace of the received signal with the three returns 
separated in time. The initial, large amplitude (clipped 
in picture) signal is the direct wave which is received 
immediately. The reflected signal due to normal incidence 
of the side lobe appears at approximately 820 microseconds, 


and the backscattered return from the crest arrives at 
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Begune 15. Typical oscilloscope trace of direct, 
reflected, and diffracted signals 
from smooth wedge 
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approximately 1.5 milliseconds. The square pulse above the 
received signal shows the sampling aperture for gating the 
return to be analyzed. 

The geometry of Figure 4 and the picture in Figure 15 
show how the source-receiver position limitations are 
imposed. If the source-receliver is too close to the crest, 
the backscatter will be masked by the stronger direct and 
reflected returns. If the range is kept constant, but the 
angle is increased, the reflected signal is delayed in time 
and overshadows the backscatter again. If the range is in- 
creased to more than half the length of the side of the wedge, 
diffraction from the back edge interferes with the diffraction 
of interest from the crest. To determine the position limi- 
tations numerically, data collection was conducted at angles 
ceo tO’, 20°, 30° at both 25 cm. and 30 cm. Figures 16 
through 23 show the normalized, backscattered intensity 
plotted versus frequency. It is noted that the backscatter 
for a wedge can be defined by the quantity of backscattered 


Signal strength (BSS), 


I 2 
a BS 
BSS, = 10 log[ 7 


Sie 


=] 
0 


where Ty 1s the intensity at one meter. However, Since it 


is the intensity ratio which is calculated in the signal 
processing, this is the quantity used in the comparison with 
theory. In theory, In6 1s inversely proportional to frequency 
for a wedge at high frequencies (Ref. 7). 


Sas 


0 Gap Kzo0uj ‘(v) e3ep TequUoUTrsedxe 
£9 = “9 = 6 ‘WO GZ = “1 = 21 ‘AaR}BOSYORG aAPpaemM YyROOoUS 


ZH ‘bau 


J. o€ 02 
Se ee 


“gT oanbtyg 


Ol 


OV - 


52 





-(-) AzoayA ‘(V) Beep TeQUOUTIAadxe 
- OT = 04 =e eG Ce eer ok ‘qoqqeosyoerq ebpem yRoousg “LT eanbtg 


ZH ‘bess 
OO! OL OS OE OZ Ol 
OP - 
Vv 
Vv 
v OE- 
vy a Vv sd Vv 

VV y yy ae 

Vv Oo 
a 

_ 4 V gs 
QM. 
Vv OD 

OL- 


Se 





a 


=) 
502 = 04 =O uo eG Ce O, = 2 ‘ZaqAeosyorq obpaem yRoouS 
ZH} ‘baa 
OO! OL OS O€ OZ 
t-—————_-—__}---—. $K_tm —————_-——————I 
V 
V 
7 V 
Vg V 
V 7 Y Vv 
9 V V 
W/L 


Azooy4 ‘(Vv) e3ep TeqUeWUTAadxs 
“8ST oeanbty 


Ot 


Ov - 


54 





-(-) AzZooyR ’(V) eIep Te j}USeUTAsdxs 


fS0O€ = 05 = 9 ‘wo Gz = OT ‘zaqjeosyorq ehpem yRooUS “°6T eANnhbTyT 
ZH ‘baw 
OOL OL OS O€ OZ Ol 
Vv 
V 
OV - 
Vv 
Vv 
V 
V Vv V O&s 
V Vv 
b OF 
V 
Vv O¢- 
V A. 
Y oo 
Vv 
OL- 


35 





-(-) AzosyA 4(M1) uotsusRzxs 
UTM obpam ATOZ eRep ‘4(V) eRzep Tej,UOUTIedxse 


£50 = 0, =—— 0 Wome = 0, = I ‘ZaqQReOSyoRq shbpom yRooUS 
ZH ‘bea 
Oo! OL OS O€ 
pp rt] 
V 
V 
V 
Z 0 0 oO 
oO V 
V 
vv Vv 
y Vv 


“oz eanbtg 





56 








-(-) Azoay3 ‘( 0) uotsuaeyxa 


UAITM Obpem AOF eRep ‘(V) e}ep Te }USUTASedxe 


‘ ot = 97 = 9 ‘WO O€ = et ‘ZaqqQeosyoeq ahpam yRyoouS “*TZ eAnhbty 
ZH ‘basa 
ool OL OS O€ 02 O| 
[{-——______—-—__________++ 2 
y Vv ¥ OV - 
Oo 
y V y 
v g O 8 
5 Vv = 
oO - OE 
Y : 
V Vv oth: 
v Fes 
Os 
y Oo. 
CO 
OL- 


ey, 





; ae Kz0auy (vy) ezep Tejuoutarzsedxa 
£002 = "09 = 69 ‘WO OF = “AT =I ‘AaQRZeOSYyDeq ebpem yRooUS *7Z oanbt 4 


4} 


t 


O¢- 


€P 


OL- 


58 





*(-) AzroayuA ‘(Vv) eIep TeqQuUSoWTASdxa 


{,0€ = 6 = 6 ‘WO OE = OL, = g 4‘Taqqeosyoereq ahpam yyooUS °¢z aaANnhtTy 


ZH ‘baa 


Ook OL OS O€ 0Z 
a 


Ot 


Ov - 


Oks 


ee 
I 


O 
sq 


OG 


gP 


OL- 


og 





The received signals were time averaged 1000 times in a 
200 microsecond sampling aperture to obtain a sSignal-to- 
noise ratio (S/N) ranging from 5 to 30 dB throughout the 
frequency spectrum; Figure 24 shows the relative S/N for the 
three surfaces. The noise was measured with both source 
and microphone on but with a data window corresponding to 
a region of nonscatter in the anechoic chamber. 

The results of Figs. 16-23 indicate a major consistent 
deviation from theory for the frequency range above 50 kHz. 
This can be explained by reverting to the beam pattern dia- 
grams. At the ranges of the experiment, the 200 microsecond 
processing window for the returned signal subtends a half~ 
angle of 28°. This angle of insonified area should be within 
the main lobe of the source-receiver pattern to insure symmetry 
and approximate a point source. The plot for 50 kHz approaches 
this critical condition where the rapid fall-off from the 
main lobe is within the angle of collected data. However, 
mmemeolar plots for the higher frequencies show increasingly 
asymmetrical patterns near the acoustic axis and lack of a 
well-defined lobe pattern expected at those frequencies. This 
non-uniformity and deviation from a point source lobe renders 
the data above 50 kHz unreliable. 

At the 25 cm. range, the data up to 50 kHz shows close 
agreement with theory for the low angles of 0° and 10°. 
However, as the angle 1s increased, the interference from 


the reflected signal causes an increasing oscillation of the 
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Spec avout the theoretical curve. The plots for the 30 cm. 
range show similar results, although the lower angles do 
not compare as well as at 25 cm. This can be attributed 

to the concurrent diffraction signal from the back edge of 
the plate. A plywood extension was added to the plate to 
see if the data at 30 cm. could be improved. The extension, 
with resulting reduction of back edge diffraction, improved 
the data slightly, except for one or two points, which may 
be due to the impedance mismatch between the plywood and 
aluminum. From these results, the geometry was restricted 
to frequencies no greater than 50 kHz, grazing angles of 0° 


moa , and a range of 25 cm. 


C. ROUGH SURFACE 

The backscatter from a randomly rough surface due to a 
transmitted signal at normal incidence was measured in order 
to compare to experimental and theoretical scattering results 
reported by Clay and Medwin [2]. Previous studies deal with 
surfaces with small slope roughnesses and conclude that the 
backscattered intensity will be the sum of a coherent com- 
ponent and an incoherent component with the coherent dominant 
at g << 1 and incoherent at g >> 1. The measured statistic 
of rms height, co, should be comparable to the value calculated 
from acoustical normal incidence backscatter in the range of 
g < 1 as described in Section III.B. The lead shot surface 
presented a backscattered intensity (not normalized) at nor- 


mal incidence that was almost identical in shape and deviated 
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less than 5 dB from the backscattered intensity from a smooth 
surface (Figure 25). The fact that for most frequencies the 
return from the wedge covered with lead shot yielded a slightly 
greater signal than from the smooth surface may have been 

the result of incoherent contribution of backscatter from 
multiple large slope off-axis scatterers. The similarity of 
these curves and the small dimensions of the lead shot compared 
to the sound wavelength led to the tentative conclusion that 
this surface appeared nearly "acoustically smooth" for the 
frequency range of interest. The term "acoustically smooth" 

1s used in the classical sense, g << l. 

The experiment now shifted to data collection from sur- 
faces coated with gravel to present a greater degree of 
roughness. 

Two methods were used to derive the statistical parameters 
of the surface. In the acoustical technique, the normal inci- 
dence mean-squared backscattered pressure (normalized to 
Squared pressure received from smooth surface at the same 
range) was plotted as in Figure 26. The backscattered pressure 
was measured at four regions of the rough surface, to obtain 
the mean squared pressure. Comparison of the curve with 
theory (eqs. 27, 29) was used to extract the rms height and 
correlation length. Alternatively, a small sample plate of 
gravel roughness was constructed. The heights ina 5 cm. square 
area were sampled every 0.1 mm. to provide a 2500 point data 


base for calculating the same statistical parameters. 
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Theoretically, for g << 1, the log of the relative inten- 
Sity is proportional to the roughness parameter, g. The 
roughness parameter (g) is related (see eq. 22) to the rms 
height (co), which can now be obtained. The latter parameter 
and the asymptotic value that the curve approaches enable 
calculation of the correlation length (L). The backscatter 
data were plotted against £* and estimates from data in Figure 


26 were calculated as follows, 


== oes 
0 
Tas 
nO LOG Ta. = 10 1G LOG “e 
0 
_ 2 2 
g = kos (cos 84 + cos 94) 
: 22 
for 90° backscatter, g = 4k a 
a Ze 
10 log = = 4,343 (24) 0% eee can/ 
0 e 
I 
A{10 log =) 
“0 -7, 2 
5 = slope = (-5.762 x 10 )a 
MET 
ae mceole liom onc 
g = 8.20% 105° em = 0.82 mn 
acoustical 5 
Oo = 2.23 mm 
msrmt 
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and 


AL“ 
Ge i Ze 
a (ASYMPTOTE) 20°r*, */ 
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PoP rOtTE = 0.38 


eae) 1101 sem 


rh 
il 


Zo) Cm 


Pe= 579. om” from window opening 


@ 


Pe oustical = 7.73 x 10 cm = 0.77 mm using Chas a cepssis 
tA coustical = 2.0 mm using °mMsxrmt 
Lsrmt = 4.0 mm. 


The statistical measurements performed by Jordan [11] on the 
sample gravel plate yielded values of o = 2.23 mm. and L = 4.0 
mm. 

Although the shape of the plotted data seems promising, 
the differences in the two statistical parmaeters are evident 
by the plot. The slope shown leads to o = 0.83 mm. vice surface 
measured o = 2.23 mm. Similarly, the acoustically determined 
correlation length is 0.77 mm. vice the surface measurements 
4.0 mm. The large discrepancies between measured and acous- 


tically determined statistical parameters can be 


oF 





explained by the inappropriate theoretical description of 

the rough surface and the inability to use a large number of 
scattering surfaces for the experimental determination of 
<p>*/p_*. The experiment used gravel that presented a jagged 
rough surface with steep slopes and cusp-like features, 
whereas the acoustic scattering theory from which o and L were 
calculated assume a gently changing slope for the rough 
Ssuerace. 

The experiment next proceeded to identify the effects of 
both a wedge and the rough surface on the backscattered sound. 
The backscattered signal strength (BSS .) for the rough surface 
alone 1S expressed by the equation for incoherent point 


scatters, 


Tas oR? | 
T 5 
0 RoW A 





BSS. = ie roge 


This definition of BSS assumes that the intensity diverges 

as = from source to scatterer, and rR iaon HOolnt Scatterer 
to receiver for a net divergence of Ree The point scattering 
assumption was verified for frequencies above 50 kHz at nor- 
mal incidence backscatter. This supports the assumption of 
point scattering at all lower frequencies, and at grazing as 
well as normal incidence. The total signal returned froma 


rough surface and wedge is the sum of these scattering effects, 


IRs eo Ins Ro 
I 2 I 2 
0 ro r 0 Ro A 
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The experimental data is a total intensity ratio from the 
composite surface, where the proportion due to the individual 
effects is unknown. The geometry of 0° and 10° at 25 cm. 
range was used so that the results from the wedge as reported 
in Section B could be separated from the plots with confi- 
dence that the remaining information was a result of the 
rough surface. The lack of any specific shape of the data, 
Figure 27, that might possibly equat€ to a simple frequency 
dependence, is found also in the literature [3,12]; as the 
backscatter due to sound at low grazing angles is determined 
by complex processes of multiple scatter and shadowing which 
have not been adequately described by any theory. The ampli- 
tude of this near grazing backscatter in the plots is signifi- 
cant in that it completely masks any diffraction effect of 
the wedge, which is 5-15 dB less under the same experimental 
Sonar lOns. 

This absence of identifiable wedge backscatter now re- 
gGuired the experiment to use a model which would reduce the 
rough surface backscatter relative to wedge diffraction. 
The previous conclusion of the "smooth" appearance of the wedge 
with the lead shot covering needed to be reevaluated. The 
slight roughness of that surface at low grazing angles might 
produce a combined backscatter with more equitable proportion 
between the magnitude of the wedge and the rough surface 
eter Tl. 

The low grazing, backscatterd signal from the lead shot 


plate was sampled over different time apertures to see the 
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varying responses of rough surface backscatter or the emer- 
gence of the diffracted signal from the crest of the wedge. 
The least time path (LTP) is the time for sound to travel 
directly to the crest of the wedge and return. A sampling 
window prior to the LTP analyzes scattering from the rough 
surface only; whereas a window that includes the LTP receives 
the combined scattering effects from the rough surface and 
the wedge. 

An additional data collection run obtained a greater 
frequency resolution by changing the transmitted frequency 
of the pulse and the number of samples taken per sampling 
window. Figures 28 and 29 show a fairly close agreement 
between the slope of the data and that theoretical slope 
for the wedge, backscatter. As the grazing angle becomes 
lower (0° situation, Figure 29), the data tends to diverge 
from the theoretical line as the magnitude increases from 
the rough surface backscatter. In the plot of increased 
frequency resolution (Figure 30), the sampled aperture of 
rough surface and wedge backscatter (@) follows a pattern 
of linearly decreasing magnitude with increasing frequency, 
yet the slope is much greater than the frequency relationship 
for wedge diffraction alone, which is theoretically eel é 
[7]. In Figures 32 and 33, the backscatter from the rough 
surface alone presents consistent patterns; although once 
again, no specific frequency dependence relationship can be 


modeled from the data. The Significant feature of this 
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backscatter effect is the magnitude of the rough surface 
backscatter. From both the fact that magnitude increases as 
grazing angle decreases and the relatively greater magnitude 
of the backscatter as compared to randomly rough surface 
backscatter in the literature [3,12], the results suggest 
the boundary wave effects proposed by Tolstoy [6]. Tolstoy 
describes a boundary wave generated from a rough surface 

for near grazing incidence and kh < 1 (h = separation of 
roughness elements). The magnitude of backscatter for the 
experiment suggests the possibility of a large backscattered 
boundary wave for surfaces of low height roughness, da/i, 

but large slope of the point scattering elements, o/L, for 


sound at low grazing angles. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


Experimental analysis of impulse wave backscatter has 
been accomplished in the laboratory for various physical 
models. 

The detected return from the shallow change of slope for 
a plane wedge shows close agreement with the theory proposed 
Eyes1ot and Tolstoy. 

The difference between the actual statistical parameters 
of the rough surface and those estimated from the spectrum 
of scattered energy at normal incidence is attributed to 
the jagged, cusp-like appearance of the surface; that is 
not accounted for by the theoretical description of the 
SuiG@race. 

The inability to positively separate the individual back- 
scatter identities from a composite surface is caused by the 
disproportion in magnitude between the two effects. The 
total magnitude in the combined spectrum suggests a large 
backscattered boundary wave for low grazing angle propaga- 
tion over interacting steep sloped rough surface elements. 

Additional experimental research in this topic could 
include changes or improvements in the areas of source/ 
receiver selection, physical model design, and theoretical 
description of rough surfaces. 

The ability to find a smaller source/receiver and yet 


Maintain sufficient, detectable output signal and receiver 


no 


a 





sensitivity could extend the agreement with theory above 

50 kHz. Also, the design of an electronic switching circuit 
to allow the use of a single microphone for both source 

and receiver would expand the flexibility in geometry. 

A rough surface wedge with a greater wedge slope would 
increase the backscatter effects from the wedge so that they 
are not masked by the return from the rough surface. 

Further iteration of the mathematical description of 
the rough surface might more realistically model this natural 


random surface to the experimental data. 
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